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Abstract 

The quartic force field of SO3 was computed fully ab initio using coupled 
cluster (CCSD(T)) methods and basis sets of up to spdfgh quality The effect 
of inner-shell correlation was taken into account. The addition of tight d func- 
tions is found to be essential for accurate geometries and harmonic frequen- 
cies. The equilibrium geometry and vibrational fundamentals are reproduced 
to within 0.0003 A and (on average) 1.15 cm 1 , respectively. We recom- 
mend the following revised values for the harmonic frequencies: u>i=1082.7, 
^2=502.6, o;3=1415.4, 6^4=534. cm" 1 . In addition, we have shown that the 
addition of inner polarization functions to second-row elements is highly desir- 
able even with more approximate methods like B3LYP, and greatly improves 
the quality of computed geometries and harmonic frequencies of second-row 
compounds at negligible extra computational cost. For larger such molecules, 
the B3LYP/VTZ+1 level of theory should be a very good compromise between 
accuracy and computational cost. 



*In memory of my colleague Dr. Jacqueline Libman OBM (1941-1997) 
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I. INTRODUCTION 



The importance of the sulfuric anhydride (SO3) molecule in atmospheric and industrial 
chemistry requires no further elaboration. 

Experimental studies of its vibrational spectroscopy have to contend with a number of 
problems such as the tendency of SO3 to polymerize, its very hygroscopic character, and its 
easy decomposition in the gas phase to SO2. Possibly this explains why only a relatively 
small body of rotation- vibration spectroscopy data is available. 

The only available anharmonic force field for this prototype planar XY 3 molecule is the 
work of Dorney, Hoy, and Mills (DHM), who proposed two model force fields. Not enough 
data are available for an experimental force field refinement. 

In recent years, however, the methodology for computing anharmonic force fields fully ab 
initio has developed to the point where the computed force fields are at least P,[3j suitable 
starting material for further experimental refinement and at best deliver essentially 

spectroscopic accuracy in their own right |7|-[10|. (By analogy with Boys' concept of 'chemical 
accuracy', ±1 kcal/mol, for energetics, we will arbitrarily define 'spectroscopic accuracy' here 
as ±1 cm' 1 on vibrational transition energies.) 



Recently, the present author published a calibration study |TT| on the anharmonic force 
field of sulfur dioxide. It was found there that the computed geometry and harmonic frequen- 
cies were critically dependent on the presence of high-exponent 'inner polarization functions' 
in the basis set, whose contribution is actually much more important than that of inner-shell 
correlation. (They also contribute as much as 8 kcal/mol to the total atomization energy 
of SO2 |rTi|12[ . This is an extreme version of a phenomenon that appears to occur more 
generally in second- row compounds [|T3| , p!4| .) Our best computed force field agreed to within 
0.0004 A and 0.03 de grees with experiment for the geometry: the errors in the fundamentals 
of SO2 were +3.9, -0.4, and +0.4 cm" 1 . It would therefore appear that the same level of 
theory would be sufficient to produce a reliable force field for SO3 as well. 

Some other recent ab initio anharmonic force field calculations on planar XY 3 systems 
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include the work of Botschwina and coworkers []15[ on CF3 and CF3 , of Schwenke |T6| on 
CH 3 , of Green et al. [17j on SiFLt, of Martin and Lee Jl8| on BH 3 , and of Pak and Woods 
|TJ on BF 3 and CF^. 



Previous ab initio studies of the vibrational force field were limited to the harmonic 
part (e.g. Pm , which used scaled quantum mechanical (SQM) |2T| techniques as well as 
MCSCF calculations). To the author's knowledge, this paper presents the first accurate 
anharmonic force field for SO3 obtained by any method. It will also be shown that the 
harmonic frequencies derived from experiment by DHM need to be substantially revised. 

II. METHODS 



All electronic structure calculations were carried out using MOLPRO 96.4 [22] with the 



tripu and scfprO patches |23| installed, or a prerelease version of MOLPRO 97.3 p4 
running on a DEC Alpha 500/500 workstation and an SGI Origin 2000 minisupercomputer 
at the Weizmann Institute of Science. 



As in our previous work on SO2 |TI[, the CCSD(T) electron correlation method p5p6 



as implemented by Hampel et al. ^7|, has been used throughout. The acronym stands for 



coupled cluster with all single and double substitutions [28] augmented by a quasiperturba- 



tive account for triple excitations ||25|| . From extensive studies (see [^9| for a review) this 
method is known to yield correlation energies very close to the exact n-particle solution 
within the given basis set as long as the Hartree-Fock determinant is a reasonably good 
zero-order reference wave function. From the value of the T\ diagnostic |30| calculated for 
SO3, 0.018, we see that this condition is definitely fulfilled. 

Calculations including only valence correlation were carried out using correlation con- 
sistent polarized n-tuple zeta (cc-pVnZ) basis sets of Dunning and coworkers |31| , |32[| . For 
sulfur /oxygen, the cc-pVTZ and cc-pVQZ basis sets correspond to (15s9p2<il//10s5p2(il/) 
and (I6sllp3d2flg/I2s6p3d2flg) primitive sets contracted to [5s4p2dlf /4s3p2dlf] and 
[Qs5p3d2flg /5sAp3d2flg], respectively. Because of the strongly polar character of the SO 
bonds, we have also considered the aug-cc-pVnZ basis sets p3|j34][, which consist of cc-pVnZ 



basis sets with one low-exponent 'anion' function added to each angular momentum. In the 
interest of brevity, the standard acronyms cc-pVnZ and aug'-cc-pVnZ will be replaced by 
VnZ and AVnZ, respectively. 

Furthermore, we considered the addition of a tight d function with exponent ct£)=3.203 
13] to the sulfur basis set; its presence is indicated by the notation VTZ+1, VQZ+1, and 



the like. 

As in our previous studies on SO and SO2 flTT] , H 2 SiO [Pj , and various second-row 



diatomics p5j[T3l , core correlation was included using the Martin- Taylor basis set |p6j^7 



This is generated by completely decontracting the cc-pVTZ basis set and adding a single high- 
exponent p function and two even-tempered series, one of three tight d functions and another 
of two tight / functions. The exponents are defined as 3 n a (rounded to the nearest integer 
or half-integer), where a represents the highest exponent of that angular momentum already 
present in the cc-pVTZ basis. Obviously, such a basis set already amply covers the 'inner 
polarization' region as well. For brevity, CCSD(T) /Martin- Taylor calculations with only the 
valence electrons correlated are denoted CCSD(T)/MTnocore, while in CCSD(T)/MTcore 
calculations, all orbitals except for the sulfur (Is) -like orbital (which lies too deep to interact 
appreciably with the valence shell) have been correlated. 

Geometry optimizations were carried out by repeated multivariate parabolic interpolation 
with a step size of 0.001 bohr or radian, and a convergence threshold of about 10 -5 bohr 
or radian. Quartic force fields were set up by finite differentiation in symmetry-adapted 
coordinates 

(1) 
(2) 
(3) 
(4) 
(5) 
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in which the are the SO bond distances, the d{ are the OSO bond angles opposite the 
corresponding r^, and S represents the out-of-plane motion. In order to keep higher-order 
contamination in the quartic portion of the force field to a minimum, fairly small step sizes 
0.01 A and radian were used and CCSD(T) energies converged to essentially machine pre- 
cision. Generation of the displaced Cartesian geometries and transformation of the internal 
coordinate force field to Cartesian coordinates were carried out with the help of the INTDER 
| 38fl program. The anharmonic spectroscopic analysis was carried out by standard second- 



order rovibrational perturbation theory using a modified version of SPECTRO p9|j40l . No 
Fermi resonances needed to be accounted for, but some of the rotation-vibration interaction 
constants needed to be deperturbed for a Coriolis resonance u> 2 ~ ^ 4 around the b axis. 

III. RESULTS AND DISCUSSION 

Computed and experimentally derived bond distances and harmonic frequencies are given 
in Table |, while computed and observed fundamentals can be found in Table ||. 

High- resolution values are available for all fundamentals except u±, for which the Bondy- 
bey and English (BE) |^T[ value of 1068.6 cm -1 appears to be the most reliable one available. 
^3=1391.5205 cm -1 was taken from the work of Henfrey and Thrush (HT) ||42|| . The rota- 
tional fine structures of the v 2 and bands overlap as well as exhibit Coriolis resonance: 
Kaldor et al. (KMDM) resolved this spectrum to give ^2=497.55 and z/4=530.18 cm -1 . 
They also definitively refuted a suggestion by Thomas and Thompson (TT) [Q that the 
assignments be reversed: the TT band origins (with the correct assignment) are z/ 2 =498.5 
and z/ 4 =529.16 cm -1 . A high-resolution study by Ortigoso, Escribano, and Maki (OEM) 
H finally yielded u 2 =497. 5679(1) and z/ 4 =530.0863(l) cm" 1 . In the process, OEM also ob- 



tained improved values of the ground state rotational constants, and particularly a revised 
r e =l. 41732 A, which is near the lower limit of the earlier DHM value, r e =1.4184±0.0010 A. 

Comparing the observed with the calculated fundamentals, it is immediately seen that 
the CCSD(T)/VTZ fundamentals leave a lot to be desired: errors with respect to the ex- 
perimental values (BE,OEM,HT) are -25.1, -15.5, -20.4, and -16.1 cm -1 , respectively, for 



vi, v ii v Zi an d ^4- Switching to a VTZ+1 basis set (i.e. adding the tight d function on S) 
dramatically cuts these errors to -1.6, +1.0, +4.8, and -2.0 cm -1 . 

Turning to Table |, we see that in addition r(SO) is shortened by no less than 0.0116 
A, thus cutting the discrepancy with experiment by more than half. Actually, the tight d 
function appears to be even more important for the quality of the results than increasing 
the basis set from VTZ to VQZ. Adding the inner polarization function to the VQZ basis 
set has a smaller effect than for the VTZ basis set (as it well should, since the VQZ basis 
set contains tighter d functions than its VTZ counterpart), but still affects the bond length 
by -0.0050 A. The main difference between the VTZ+1 and VQZ+1 harmonics is that uj^ 
drops by about 4.5 cm -1 , removing the largest error remaining in the fundamentals at the 
CCSD(T) /VTZ+1 level. 

Given the strongly polar character of the S-0 bonds in SO3, the addition of anion func- 
tions is expected to have a nontrivial effect. From VTZ+1 to AVTZ+1, the harmonics 
are lowered by no less than -13.2, -5.7, -20.1, and -8.0 cm -1 , respectively; from VQZ+1 to 
AVQZ+1, the effect is much weaker: -3.6, -2.4, -6.0, and -2.4 cm -1 . While the addition of 
diffuse functions lengthens r(SO) by +0.0030 A between VTZ+1 and AVTZ+1, this likewise 
becomes much less significant from VQZ+1 to AVQZ+1 (+0.0009 A). 

Computing anharmonic corrections with these large basis sets for a four-heavy atom 
molecule is beyond the presently available computational resources. However, it is seen 
here that even from CCSD(T)/VTZ to CCSD(T)/VTZ+1, the basis set effect on the anhar- 
monic corrections is very modest (-0.10, -0.27, -0.07, and +0.01 cm -1 ), while for SO2 [[11 



their basis set dependence was likewise found to be quite mild. We therefore opt for the 
CCSD(T)/VTZ+1 anharmonicities and will merely substitute an improved bond distance 
and improved harmonic frequencies in the analysis. 

Doing the latter with the CCSD(T)/AVQZ+1 data yields harmonic frequencies which 
are systematically too low, with r e of course still being 0.0053 A longer than experiment. 

A complete CCSD(T)/MTcore harmonic frequency calculation proved too demanding 
in terms of computational resources: what could be succesfully completed was a univariate 
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optimization of r e , with u\ being obtained as a by-product. The CCSD(T)/MTcore r e falls 
within the error bar of the DHM value and is just 0.00032 A longer than the OEM value. 
The contributions of inner-shell correlation to r e and u)\ turn out to be much more modest 
(—0.00283 A, +4.2 cm -1 ) than those of inner polarization. 

If we consider the ratios c<jj(MTcore)/c(j i (MTnocore) for u>i and for the harmonic fre- 
quencies of SO2, we find that these ratios are remarkably similar: 1.0039, 1.0037, 1.0038, 
and 1.0036. (They correspond almost exactly to [r e (MTcore)/r e (MTnocore)] 2 .) Under these 
circumstances, we might be able to obtain 'best estimate' harmonic frequencies by scaling 
the CCSD(T)/AVQZ+1 values by cJi(MTcore)/c;i(MTnocore) = 1.00386, which leads to the 
following values: 1081.2, 503.1, 1415.1, and 531.6 cm" 1 . 

We now substitute these harmonic frequencies and the CCSD(T)/MTcore geometry in 
the spectroscopic analysis, and thus obtain the fundamentals labeled 'best' in Table |I[ 
Compared to the (BE,OEM,HT) set of experimental values, the remaining discrepancies 
are -1.5, +0.4, -0.3, and -2.4 cm" 1 , or 1.15 cm" 1 on average. We can thus safely claim 
'spectroscopic accuracy' for our best force field. 

Turning now to the rotational constants (Table |111[ ), we see that the computed and ex- 
perimentally derived (OEM) B e agree to -0.02 %: C e brings no additional information since 
it is fixed by the planarity relation C e = B e /2. From our computed rotational constants we 
can however now determine computed B and Cq, which likewise turn out to be in excellent 
agreement with the observed values: -0.06 % and -0.03 %, respectively. (B and C are 
independent data because of the inertial defect.) The fact that the rotational constants are 
consistently computed slightly too small is consistent with our bond distance being slightly 
longer than the true value. From the relationship (in this case) tq = r e ^J ' B e j Bq, we find 
r (calc.)=l. 42004 A, compared to an experimental value ro(obs.)=l. 4198+0. 0002 A, 
and one derived from the OEM B of 1.41963 A. The discrepancy of +0.0004 A between the 
computed and OEM-derived r$ is consistent with the discrepancy of +0.0003 A between the 
computed and OEM r e values. Aside from the observation that this would be considered 
excellent agreement between theory and experiment by any reasonable standard, this sug- 
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gests that the OEM equilibrium bond distance of 1.41732 A would be accurate to 0.0001 A 
or better. The computed r g (electron diffraction) distance from our force field, 1.42275 A, 



is substantially longer than the experimental electron diffraction result | 47fl , 1.4188±0.003 
A. Since this latter study also finds an unrealistically short r e =1.4142 and r z =lA16g A(our 
own computed r z =l. 42143 A), we can safely conclude that the electron diffraction result is 
in error. (For a review of the different types of bond distances discussed, see Ref. PSfl .) 

Our best harmonic frequencies agree relatively well with those obtained by DHM from 
the experimental fundamentals and a valence model for anharmonicity, except for lu\ which 
still is about 16 cm -1 too low in the better of their two models. (The discrepancies for 
co>3 and cl>4 are still substantially bigger than the differences between the computed and 
observed fundamentals, and it can safely be stated that the present ab initio values are 
considerably more reliable than the experimentally derived ones.) It thus comes as no surprise 
that the anharmonicity constants (Table [IV|) obtained by DHM from their valence models 
(particularly uii— z/ 1; which has the wrong sign) differ profoundly from the presently computed 
set, which clearly is the more reliable one. 

A set of 'experimental' harmonic frequencies can be derived from our best force field and 
the observed fundamentals by iteratively substituting oo\ n+1 ^ = cj] n ' + z/j cxpt, l — yf 1 in the spec- 
troscopic analysis. The values thus obtained are given as the entry labeled 'Recommended' 
in Table |. 

Coriolis and rot at ion- vibration coupling constants can be found in Table |TV| . The com- 
puted Coriolis coupling constant -B e C^ is found as 0.1764 cm -1 from our best force field, 
which agrees reasonably, but not very, well with the OEM 'Fit I' value of 0.191694(460) 
cm^ 1 . They note that their fit exhibits very strong dependency between -B e £|| and such 
parameters as the rotational Z-doubling constant q± and the rovibrational coupling constants 
ct2B and a4B- Upon constraining B e C,^ to the force-field derived value of 0.1801, they found 
(Fit II in their paper) that all the abovementioned constants change drastically. It is note- 
worthy that our computed a 2 B=0. 00031, a 4 £=-0. 00050, and g 4 =0. 00063 cm -1 agree much 
better with the 'Fit II' values of 0.000150, -0.000404, and 0.000497 cm" 1 , respectively, than 
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with the 'Fit F values of -0.000381, -0.000140, and -0.000047 cm" 1 , respectively. Linear ex- 
trapolation suggests that much of the residual discrepancies between our computed and the 
'Fit IF would disappear if our value had been substituted in the experimental analysis. 

OEM also obtained 0.2c and a^c values, which are in excellent agreement with our calcu- 
lations, as are the a 3 B and a 3 c values of HT. Our computed a\s and ol\c are larger than the 
model-derived values of DHM, which appear to be on the low side for all other constants as 
well. The centrifugal distortion constants are quite small, and agree with the experimental 
values of OEM to within the latter's uncertainties. 

The trends in the computed quadratic force constants (Table |V|) closely parallel those 
in the harmonic frequencies. (The most striking difference with the DHM quadratic force 
constants, which reproduce the fundamentals rather than either of their sets of harmonic 
frequencies, lies in the stretch-bend coupling constant F 34 .) 

In order to stimulate further research on the vibrational spectrum of SO3, the symmetry- 
unique cubic and quartic force constants in symmetry coordinates have been made available 
in Table |VT|. The force fields in Cartesian, symmetry, and normal coordinates can also be 
downloaded in machine-readable form on the World Wide Web at the Uniform Resource 
Locator (URL) http : //theochem. weizmann. ac . il/web/Papers/ so3 .html 



IV. A NOTE ON LOWER-LEVEL CALCULATIONS 

Some readers might wonder how well less computationally intensive methods would do 
for the mechanical properties of SO3 and the other sulfur oxides, and whether the presence 
of the tight d functions is still relevant at that accuracy level. 

In order to answer these questions, we have carried out geometry optimizations and 
harmonic frequency calculations for SO, S0 2 , and S0 3 using the popular B3LYP density 
functional method as implemented in GAUSSIAN 94 |5g. The VTZ, VTZ+1, and 



AVTZ+1 basis sets were considered, as was the popular 6-31+G* basis set. The results are 
summarized in Table |V11[ 

It is immediately seen that the 6-31+G* basis set systematically overestimates bond 
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lengths by no less than 0.035 A, and (largely as a result thereof) underestimates stretching 
frequencies by as much as 80 cm -1 and the SO3 out-of-plane bending frequency by about 50 
cm -1 . These errors are substantially reduced by using the VTZ basis set. However, at very 
small additional expense, the addition of a tight d function on S leads to quite respectable 
agreement with experiment: residual errors for the B3LYP/VTZ+1 harmonic frequencies are 
+5.7 cm" 1 in SO, {+0.4,+16.2,+6.5} cm" 1 in S0 2 , and {-5.7,-8.4,-10.5,-6.4} cm" 1 in S0 3 . 
Trends in the effect of the tight d function closely parallel those seen at the CCSD(T) level, 
which is not surprising since it is essentially an SCF rather than a dynamical correlation 
effect. |TJjrj 

At the B3LYP/VTZ+1 level, all bond lengths are now within +0.006-0.008 A of experi- 
ment, and the OSO angle is in excellent agreement with experiment. Overall, performance 
with the VTZ+1 basis set is as good as we can reasonably hope to get (e.g. ||52|| ) at the 
B3LYP level. 

It is therefore clear that the addition of tight d functions to the basis set is eminently 
worthwhile even for less than 'benchmark quality' calculations on second-row compounds, as 
well as that B3LYP/VTZ+1 would represent an excellent compromise between accuracy and 
computational cost for geometry and frequency calculations on larger second-row systems. 

V. CONCLUSIONS 

The first-ever accurate anharmonic force field for SO3 has been obtained fully ab initio. 
We have been able to establish that: 

• r e is reproduced to within +0.0003 A, and the fundamentals to within 1.15 cm -1 , on 
average; 

• like for SO2 and (to a lesser extent) for second-row compounds in general, the addition 
of tight d functions ('inner polarization functions') to the basis set is essential for 
accurate results; 
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• the following revised values are recommended values for the harmonic frequencies: 
cji=1082.7, ^=502.6, w 3 =1415.4, u 4 =53A.O cm" 1 ; 

• our computed rovibrational coupling, rotational /-doubling, and Coriolis coupling con- 
stants suggest a preference for the set of constants in 'Fit II' in Ref. JR| (OEM) over 
those in 'Fit I'. 

In addition, we have shown that the addition of inner polarization functions to second- 
row elements is highly desirable even with more approximate methods like B3LYP, and 
greatly improves the quality of computed geometries and harmonic frequencies of second- 
row compounds at negligible extra computational cost. For larger such molecules, the 
B3LYP/VTZ+1 level of theory should be a very good compromise between expense and 
accuracy. 
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TABLES 



TABLE I. Equilibrium bond distance (A) and harmonic frequencies ( 


3H1 ) 


of SO 


3 






r e 


u\(a' x ) 




4) 






e') 


CCSD(T)/VTZ 


1.43753 


1057.7 


487 


5 


1395 


1 


518 


1 


CCSD(T)/VTZ+1 


1.42594 


1081.1 


503 


6 


1420 


2 


532 





CCSD(T)/VQZ 


1.42780 


1071.1 


496 


2 


1405 


6 


526 


8 


CCSD(T)/VQZ+1 


1.42279 


1080.6 


503 


6 


1415 


7 


532 





CCSD(T)/AVTZ 


1.44038 


1043.9 


482 


4 


1374 


1 


509 


7 


CCSD(T)/AVTZ+1 


1.42890 


1067.9 


497 


9 


1400 


1 


594 


n 

u 


CCSD(T)/AVQZ+1 


1.42372 


1077.0 


501 


2 


1409 


7 


529 


6 


CCSD(T)/MTcore 


1.41764 


1092.5 














CCSD(T)/MTnocore 


1.42047 


1088.3 














Best estimate" 


1.41764 


1081.2 


503 


1 


1415 


1 


531 


6 


Experiment 


1.41732 |4E|,1.418 4 [| 
















£ VF' model (DHM) Q 




1048.08 


503 


81 


1408 


96 


538 


64 


'Extended' model (DHM) Q 




1064.89 


505 


97 


1410 


00 


535 


62 


Recommended 6 




1082.7 


502 


6 


1415 


4 


534 






(a) harmonics obtained by scaling CCSD(T)/AVQZ+1 values with ratio of 1.00386 between 
CCSD(T)/MTcore and CCSD(T)/MTnocore value for u) 1 (see text) 

(b) this work. Obtained by iteration of harmonics with best computed force field to exactly 
reproduce experimental fundamentals. 
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TABLE II. Anharmonic corrections ( 


cm 1 ) and fundamentals (cm 1 ) 


for S0 3 








UJ 3 - Z^ 3 


W4 — ^4 


CCSD(T)/VTZ 


14.182 


5.342 


23.989 


3.911 


CCSD(T)/VTZ+1 


14.083 


5.067 


23.921 


3.920 


Best estimate" 


14.143 


5.069 


23.930 


3.921 


'VF' model (DHM) 


-19.9 


6.3 


19.1 


8.5 


'Extended' model(DHM) 


-3.1 


8.4 


20.1 


5.4 




v\ 


^2 


^3 




VTZ 


1043.5 


482.1 


1371.1 


514.0 


VTZ+1 


1067.0 


498.6 


1396.3 


528.1 


Best calc. a 


1067.1 


498.0 


1391.2 


527.7 


Expt. 


1068.6 6 


497.5679(l) c 


1391.5205** 


530.0863(l) c 


(a) from substituting CCSD(T)/MTcore geometry and 'best estimate' harmonic frequencies 
in spectroscopic analysis of CCSD(T)/VTZ+1 force field 



(b) Ref. (BE) 

(c) Ref. H (OEM) 

(d) Ref. H (HT) 
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TABLE III. Computed and observed bond distances (A), rotational constants (cm 1 ), and 
centrifugal distortion constants (cm -1 ) of SO3 





Best calc. 


Expt. 




1.41764 


1.41732 a , 1.4184±0.0010 fe , 1.414 2 c 


ro 


1.42004 


1.4198±0.0002 d , 1.41963 6 


r 9 


1.42275 


1.4188±0.003 c 


r z 


1.42143 


1.416 9 c 


B e 


0.34962 


0.34968 a , 0.34923 fc 


C e 


0.17481 


Be/2 


B 


0.34844 


0.3485439 a , 0.34857 6 


Co 


0.17393 


0.173984 a , 0.17402 6 


W 7 Dj 


3.092 


3.096(8) a 


W 7 D JK 


-5.452 


-5.47(2) a 


W 7 D K 


2.543 


2.55-f 



(a) Ref. m (OEM) 

(b) Ref. § (DHM) 

(c) Ref. 

(d) Kaldor and Maki @ 

(e) From B of OEM 

(f) from planarity relation Dk 



-{2D j + 3D JK )/A 
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TABLE IV. Anharmonicity constants, rotation-vibration coupling constants, and rotational 



/-doubling 


; constants of SO3. 


All values are in cm 1 . 


Constants marked with 


an asterisk have been 


deperturbed for Coriolis resonance 










'VP 1 ' mndpl 


'Extended' model 


Experiment 




/ f" rl 1 nmfV 1 
I L 11 la VVUIri 1 


DHM 


DHM 


OEM 


&1B 


n ooosi 


00067 


0.00067 




a 2B 


U.UU4UO 




0.00425 




a 2B 


n nnn^i 






-0.000381 a , 0.000150 6 


a 3B 


n nm 1 1 


o nm o^ 

U.UU1UO 


0.00103 


0.001132(l) c 


&4B 


-U.UUZDo 


fl 00971 
-U.UUZ I 1 


-0.00269 




a 4B 


-f) 00050 






-0.000140 a , -0.000404 6 


a lC 


n nnnAi 


000^*3 


0.00033 




OL2C 


n nnm ^ 
-u.uuuio 


00096 


-0.00019 


-0.000130^ 


a 3C 


ono^k 


000^4 


0.00054 


0.0005999(3) c 


a 4C 


0001 5 


0000^ 


0.00008 


0.000157 d 


13 




00009 


0.00002 




14 


n nn^nn 

U.UUOUU 


00*34^ 


0.00345 




ll 


n nnofi^ 






-0.000047 a , 0.000483 6 


-°S24 


0.1764 






0.1917(5) a , 0.1801 6 


^11 


9 n^i 


-l.oi 


-1.38 




X\2 


1 113 


in Sfi 


5.08 




X\3 


-9.388 


-7.02 


-6.94 




X , , 


-1 249 


24 27 


10.27 




X22 


n 7A9. 


-4 24 


-3.21 




X23 


-4.403 


-1.22 


-2.36 




X24 


0.284 


-1.98 


-2.18 




X33 


-5.434 


-4.66 


-4.63 




X34 


-3.553 


1.00 


-1.60 




X44 


-0.011 


-6.20 


-2.63 




G33 


2.828 


2.46 


2.47 




G34 


-0.279 


-0.20 


-0.18 




G44 


0.150 


6.25 


2.08 




i?44 


2.528 









(a) Ref. @ (OEM), Fit I. 

(b) OEM, Fit II (5(24 constrained to force field value 0.1801 cm- 1 ). 

(c) Ref. H (HT) 

(d) Ref. H (OEM) 
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TABLE V. Quadratic force constants for SO3 in symmetry-adapted internal coordinates. Units 
are aJ, A, and radian, and the restricted summation convention is used. 





F n 


F22 


F33 


F34 


F44 


DHM a 


5.374 


0.309 


5.269 


-0.460 


0.617 


CCSD(T)/VTZ 


5.27167 


0.30839 


5.24224 


-0.33254 


0.61456 


CCSD(T)/VTZ+1 


5.50748 


0.32393 


5.44096 


-0.36341 


0.63755 


CCSD(T)/AVTZ 


5.13533 


0.30320 


5.08914 


-0.33047 


0.59721 


CCSD(T)/AVTZ+1 


5.37410 


0.31787 


5.29261 


-0.36336 


0.62080 


CCSD(T)/VQZ 


5.40597 


0.31526 


5.32635 


-0.35288 


0.62686 


CCSD(T)/VQZ+1 


5.50201 


0.32249 


5.40691 


-0.36436 


0.63458 


CCSD(T)/AVQZ+1 


5.46613 


0.31990 


5.36281 


-0.36459 


0.62976 


CCSD(T)/MTcore 


5.62424 










CCSD(T)/MTnocore 


5.58098 











(a) reproduce older values of fundamentals, not harmonic frequencies 
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TABLE VI. Computed cubic and quartic force constants for SO3 in symmetry-adapted internal 
coordinates. Units are aJ, A, and radian, and the restricted summation convention is used. 





CCSD(T)/VTZ 


CCSD(T)/VTZ+1 


Fill 


-7.31818 


-7.62772 


^122 


-0.58764 


-0.60070 


Fl3a3a = ^13636 


-20.91209 


-21.57875 


Fl4a4a = ^14646 


-0.98244 


-1.01975 


Fl3a4a = -^13646 


1.02568 


1.05182 


i*3a3a3a = — ^3o363b/3 


-4.84233 


-4.99590 


F3a3a4a = — ^3o3i)4b/2 = — -p36364a 


-0.29699 


-0.32604 


F3a4a4a = — -F3a4646 = — ^36464a/ 2 


0.37460 


0.38143 


-F4o4o4a = — -P4a464b/3 


-0.15162 


-0.15666 


^1111 


6.05065 


6.12886 


-Fll22 


0.40348 


0.39441 


-F2222 


0.06145 


0.06940 


-Fll3a3a = -Fll3i>3(> 


33.41738 


33.72571 


F\V3a4a = ^113646 


-1.13882 


-1.10147 


-Fll4a4a = -^1 14646 


0.68428 


0.71380 


-F223a3a = -^223636 


-0.29300 


-0.30274 


-F223o4a = -^223646 


-0.34408 


-0.36252 


-F224o4a = -^224646 


0.26088 


0.27144 


-Fl3a3a3a = — -Fl3a3()3b/3 


15.54250 


15.63856 


-Fl3o3o4a = — -Pl36364a = ~ -Fl3a364!>/2 


0.70144 


0.78582 


Fl3a4a4a = — -Fl3a4f>46 = ~ -Fl364f>4a/ 2 


-0.39367 


-0.47240 


-Fl4a4a4a = — -Pl4a4646/3 


0.25712 


0.26586 


^30303030 = F3a3a3636/2 = -F3i,3f,3(,3f, 


8.41884 


8.50409 


-F3o3o3a4a = ^3a36364a = -P3a3a3646 = ^36363646 


0.04720 


0.08160 


^3o3o4o4o = ^36364646 


0.22378 


0.21243 


-F3a3a4646 = ^36364o4o 


-0.36425 


-0.40564 


-P3o364a46 = 2(F3 a 3 a 4 a 4 a — -F3 a 3 a 464i>) 






F3a4a4a4a = ^36464a4a = ^3a4a464i> = ^36464646 


-0.27446 


-0.29029 


F4a4a4a4a = F iaiaibib /2 = F ibibibib 


0.12569 


0.12956 
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TABLE VII. Basis set convergence for geometries (A, degrees) and harmonic frequencies (cm 1 ) 
of S0„ (n=l-3) at the B3LYP level 







6-31+G* 


VTZ 


VTZ+1 


AVTZ+1 


Experiment 


SO 


r e 


1.5157 


1.4998 


1.4891 


1.4888 


1.48108 H 






1118.3 


1146.6 


1157.4 


1156.2 


1150.695(8) H 


so 2 


r e 


1.4655 


1.4504 


1.4367 


1.4380 


1.43076(13) §§ 




0e 


118.68 


118.31 


119.25 


119.16 


119.33(1) §4} 






1133.0 


1164.7 


1184.1 


1177.8 


1167.91(4) @ 




u> 2 


498.4 


516.5 


522.6 


519.6 


522.21(3) [55 


1 






1315.9 


1352.1 


1388.3 


1376.2 


1381.82(2) H 


so 3 


r e 


1.4543 


1.4389 


1.4259 


1.4270 


1.41732 H 






1018.1 


1051.2 


1077.0 


1071.4 


1082. T 




LLl-2 


452.9 


476.3 


494.2 


492.2 


502.6 a 






1334.7 


1375.3 


1404.9 


1394.1 


1415.4 a 




UJ4 


492.5 


512.2 


527.6 


523.8 


534.0 a 



(a) Recommended values (this work). 
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